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Abstract 

We solve for the spectrum of the Laplacian as a Hamiltonian on R'^ — D and in 
— B. A self-adjointness analysis with 9B and dM as the boundary for the two 
cases shows that a general class of boundary conditions for which the Hamiltonian 
operator is essentially self-adjoint are of the mixed (Robin) type. With this class of 
boundary conditions we obtain "bound state" solutions for the Schroedinger equation. 
Interestingly, these solutions are all localized near the boundary. We further show 
that the number of bound states is finite and is in fact proportional to the perimeter 
or area of the removed disc or hall. We then argue that similar considerations should 
hold for static black hole backgrounds with the horizon treated as the boundary. 

1 Introduction 

Ever since Bekenstein's suggestion that black holes are thermodynamical objects that have 
entropy [T] and Hawking's discovery that black holes radiate [2], there have been numerous 
attempts to explain the thermodynamical properties of black holes from more fundamental 
quantum principles. The celebrated area law for black hole entropy suggests that the 
microscopic degrees of freedom reside on the black hole horizon and that the number 
of these states is proportional to the exponential of the horizon area. It is generally 
believed that a full understanding of black hole thermodynamics would require knowledge 
of quantum theory of gravity. Black hole entropy has been calculated in several ways - 
(i) brickwall (ii) string theory, and (iii) loop quantum gravity (LQG). It is well known in 
many of these formulations that states contributing to the entropy are obtained from a 
conformal field theory near the boundary [21 HI [5l |6l [7] . 

Although entropy calculations in string theory and LQG give results in agreement with 
what is expected from semiclassical arguments, they require details of quantum gravity. 
Nevertheless, it is also believed that, at least for large black holes, it must be possible to 
explain black hole entropy without requiring the details of quantum gravity. One such 
proposal is that the black hole entropy is the result of quantum entanglement between 
the degrees of freedom of a scalar field that are outside and those which are inside the 
event horizon (since the degrees of freedom inside the horizon are inaccessible to a distant 
observer), see [8j. 

't Hooft has proposed that black hole states form a discrete spectrum and has shown 
that the number of energy levels that a particle can occupy in the vicinity of a black hole 
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is finite if one imposes a cut-off in the form of a brick wall near tlie horizon [9l [TO] . Based 
on an earlier work [11], Bekenstein and Mukhanov [12] independently suggested that the 
black hole discrete mass levels give rise to the black hole entropy as well as the thermal 
radiation from it. These calculations, though interesting, are not illuminating in explaining 
the origin of the discrete spectrum of energy levels. 

Taking seriously the view that entropy and Hawking radiation of large black holes 
should not depend critically on the details of quantum gravity, we suggest an alternate 
proposal. The main ingredient behind the new proposal is the existence of bound quantum 
states in black hole backgrounds resulting from the self-adjoint extension of the Laplacian 
on these backgrounds. If these bound states are localized near the horizon and if their 
number is finite, then one has a candidate model to explain black hole entropy. Possibility 
of bound state solutions in the vicinity of black hole horizon has been explored previously 
also jI3j . 

Here we carry that program further but now restricting ourselves to the Schroedinger 
equation and its bound state solutions. The crucial requirement is that the boundary con- 
ditions are determined through self-adjointness. This serves two purposes: (i) the evolution 
for arbitrary combinations of modes will be unitary, and (ii) arbitrary functions in the 
domain of the Laplacian can be expanded in a complete set of eigenmodes. Interestingly, 
we find that most of the expectations mentioned in previous paragraph are realized. We 
find that there is a finite number of discrete bound state solutions and that these solutions 
are localized near the horizon. We also find that the number of these bound state solutions 
is proportional to the horizon area. 

To present the key ideas and the intuition involved in this proposal we start, in Sec. 
|2l by considering toy models in 2 -|- 1 and 3-1-1 dimensional fiat spacetimes with a ball 
removed from the spatial sections. We solve for the spectrum of the Laplacian and look 
for bound state solutions. To figure out the suitable domain of the operator in which it is 
self-adjoint we perform self-adjointness analysis. We find that the boundary conditions for 
which this operator is self-adjoint are the so-called Robin (mixed) boundary conditions. 
For the most general boundary conditions in which the operator is self-adjoint see [14j . 

We find that the Hamiltonian admits bound state solutions in the exterior of this model 
black hole and that the number of these bound state solutions is proportional to the "area" 
of the horizon (circumference of the disc in 2 + 1 dimensions). In other words, for these 
model black holes, we find that the entropy is proportional to the area. 

After this warm-up, in Sec. [3l we apply similar ideas to the case of 2 + 1 dimensional 
nonrotating Banados, Teitelboim and Zanelli BTZ black hole. However, due to the non- 
trivial metric and the existence of horizon (a null surface), the analysis now becomes much 
more subtle (see [15] for instance). After identifying the Laplacian we proceed with its self- 
adjointness analysis (using measure appropriate for this background). We show explicitly 
that the Laplacian is self-adjoint with mixed Robin type boundary conditions in tortoise 
coordinates. 

One important feature to note here is that the use of mixed boundary conditions, as 
opposed to the usual Dirichlet or Neumann boundary conditions, naturally introduces a 
length parameter in the problem. In the present context, where one is taking into account 
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both the gravitational and the quantum effects, it is natural to take this parameter to be of 
the order of the Planck length Ip. The Dirichlet and the Neumann boundary conditions turn 
out to be two special cases of the mixed boundary condition where the length parameter 
goes to infinity and zero, respectively. The effect of these boundary conditions has been 
studied in an analysis of the billiards problem in [16]. In fact these modes in a compact 
region of space give phenomena like whispering gallery modes and are also responsible for 
evanescent waves. That boundary conditions and self-adjoint extensions are in one to one 
relation was established by Asorey et. al. [H]. We will comment on this later in the 
discussion. Interesting new constraints arise when we have Dirichlet conditions at isolated 
points, see Berry et. al. [I6]. It is worth noting that the occurrence of localized bound 
states in the presence of mixed boundary conditions have implications for edge states in 
the context of quantum hall effect as well. 



2 Toy models 

We are interested in finding the bound state solutions of the time-independent Schroedinger 
equation 

- A^(r, $) = A^/'(r, $) . (1) 

Here A is the Laplacian operator, r is the radial coordinate, and $ designates all the 
angular coordinates. Note that for bound state solutions A < 0. The solution to the above 
equation is to be obtained subject to the Robin boundary condition 

Ki; + dji^ = 0, (2) 

where k is a constant and the derivative in the second term is evaluated along the outward 
normal it to the boundary. It is well known that the Laplacian operator is self-adjoint 
in the domain of functions subject to the above condition ([2]). We note that the 
constant k, has dimensions of (length)"^. This parameter would not be present if we 
use the usual Dirichlet or Neumann boundary conditions (which are the limits of Robin 
boundary condition corresponding, respectively, to k,~^ or k — )■ 0). 



2.1 2+1 dimensional fiat spacetime 

We begin by considering the case of two dimensional Euclidean plane with a disc of radius 
Tb removed from it. For two spatial dimensions ([1]), in terms of polar coordinates, becomes 

l^(r^]+l^ = X^ (3) 
r dr \ dr J 89'^ ' 

where we have made the replacement A — )■ —A so that for bound states A > 0. Since our 
boundary conditions do not mix the coordinates we can solve using the ansatz ilj{r, 6) = 
e*"^i?(r) (with n an integer) resulting in the following r dependent equation: 

d'^R IdR n^R 
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Before solving this equation we check whether the Laplacian, 

^ _ rf2 Id 

dr"^ r dr ' 

is self-adjoint [T71[l3l[T8] on the domain D(0) = {^/'|'?/' continuous and G £^[rb, oo), /t^/'(rb)- 
■?/''(rb) = 0}, defined by the boundary condition in ([2]). Here prime (') denotes derivative 
with respect to r. 

We start by confirming that the operator O is symmetric, that is, (0, Oip) — {0(j), ip) = 
V0, t/' G D{0) where the inner product, using the measure rdr, is 



(0, Oifj) = / r(j)*Oiljdr . (5) 

Substituting for O we get 

(0,O^)-(60,^) = roo(0:.^-^oo^) , (6) 

after an integration by parts and using the boundary condition in the resulting expression. 
If '000 and its derivative falloff sufficiently fast at infinity, the above boundary term is zero 
thereby confirming that the operator is symmetric. 

To check if the operator is self-adjoint we need to find the domain of the adjoint operator 
- the set of all those 0's, not necessarily belonging to D{0), which satisfy {(p.Oip) = 
{0(f), ip)^ ip E D{0). Using the same procedure as above, but with arbitrary, we get 



roo0:o ^ - f^ncpl^A, - ^000^00 ^ + niJr, ^ = • (7) 



Since and its derivative fall off sufficiently fast at infinity, the above equation gives 
following condition on at Vb. 



'^'^-= dr 



This implies that G D{0), that is, the domain of the adjoint is the same as the domain of 
the operator itself. This establishes that the operator is self-adjoint with mixed boundary 
conditions. 

Having proved the self-adjointness of the operator in (jlj), we now come to the solutions 
of that equation. These are given by modified Bessel's function /„(V^r) and Kn{y/Xr). 
Requiring the solutions to be square integrable on [rb, oo) with measure rdr {rb being the 
boundary), rules out the exponentially growing In{\/Xr) solution, and we are left with 

R{r) = cKniVXr) , (9) 

where c is a normalization constant. The complete solution is thus given by 

0;(r,^) = cir„(yAr)e^"^ (10) 
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To find the spectrum, we now impose the boundary condition in ([2]) at r = obtaining 

Here the sign on the right-hand side (rhs) is positive since the derivative is evaluated 
along the outward normal which, in the present context, points towards the origin. Since 
K'^{r) < 0, we find that there will be bound solutions only for n < 0. By plotting the left- 
hand side (Ihs) and the rhs of the above equation (for different values of n) as a function 
of a/a we find that there are only a finite number of bound solutions for a given value of 
K, and Tfe, these being given by the intersection of the above mentioned curves with the 
horizontal line representing the constant value of k, see Fig. [H 




Figure 1: 

Bound solutions in two dimensions for = 2. The horizontal line corresponds to k = —5. 

To find how the number of bound state solutions scale with the boundary radius r^, 
we start by noting that the derivative of the rhs of (ITT]) with respect to a/A in the limit 
a/a — )■ and under the assumption that n is sufficiently greater than V^, tends to zero 
(remark: it turns out that we only need n > 2 for the argument to work). This implies 
that the tangent to the curves corresponding to the rhs of (|TT|) are parallel to the a/A axis 
and therefore also to the k = constant curve. 

Thus, in the limit ^/X — )■ the curve corresponding to some order rim (> 2, see the re- 
mark above) of the modified Bessel function will either be coincident with the k = constant 
curve or will be the first curve not intersecting it at all (in either case the solution corre- 
sponding to order n = rim will not be a bound state). Thus, one can find the dependence 
of the total number of bound state solutions on the radius r?, of the removed disc. For this 
we take the limit a/A — )■ on the rhs of (ITTl) with n = rim and use the approximation 

^a(a;)~^0) iV^Tl^X>0) (12) 

to obtain 

lim v/A^^^>J^ = -^. (13) 
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Equating this to k one finds that Um scales hnearly with radius 

nm = -Krb. (14) 

From the arguments presented above, the maximum number of bound state solutions 
is equal to rim if the rhs of ( fT4l) is an integer (the total number of bound state solutions 
is Um and not rim — 1 since n = also corresponds to a bound state solution). If it is not 
an integer, then the number of bound state solutions is the smallest integer greater than 
—KTh. We also note that except for the case n = there is a two-fold degeneracy in the 
number of bound state solutions since n — )■ —n gives the same solution. 



(a) En v/s n 



20 40 
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(b) < r„ > vj s n 



Figure 2: 

Energy spectrum and the expectation value < r„ > in two spatial dimensions for = 1 
and K = —100. 

From (iTTl) it is difficult to obtain an analytic expression for the energy of bound states. 
However the spectrum can be obtained numerically and in Fig. ?? and ?? we show the 
dependence of energy En and the expectation value < r„ > on n, the order of bound states, 
for the choice k = —100 and = 1. 

From the figure we note that the energy of the lowest state Eq is close to —k'^. The 
spectrum is approximately given by En ~ —k^ + n^/r^ and thus with increasing n the 
energy approaches zero. We also note that for all n the expectation value < r„ > is always 
close to Tb showing that these states are all localized near the boundary. This is a general 
feature of the solution as can be verified by choosing different values for k, and r^,. 



2.2 3+1 dimensional flat spacetime 

We now repeat the above calculations for the case of — B^. Using the expression for the 
Laplacian in spherical coordinates, the eigenvalue problem to be solved is 



r2 Or dr)^ sin 9 09 v'"" ^ 09 ) ^ sin^ 9 0<P^ ' ^^^^ 
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where as before A > 0. The above equation can be solved easily by using separation of 
variables writing ip{r,9,(j)) = R{r)Yim{0 , (p) , where Yim{0,(j)ys are the spherical harmonics 
solving the angular part of the equation. The radial part of the equation is given by (in 
subsequent analysis we replace / by n) 

(fR , 2dR n{n + l)R 

^ + rd7 r^ = ^^- (16) 

Since the form of the Laplacian is the same as that in (jlj), it is clear that this again is a 
self-adjoint operator on the domain of interest. As in the previous case the above equation, 
with the requirement of square integrability on (r;,, oo) with measure r^c/r, is solved by the 
modified Bessel function of the second kind K^j^i {\fXr) with the only difference that now 
the order of the modified Bessel function is {n+1/2). The complete solution to the equation 
is 

V'(r, 0) = c^^^ F„^(^, 0) . (17) 

We now impose the boundary condition in ([2]) and as in the two dimensional case we 
have condition only on the radial solution 

KR{n) - R'in) = . (18) 

After substituting for R this becomes 

2k = -2VX^^ (19) 

We again find that bound solutions are possible only for k < and that for a given value 
of K and Tb only a finite number of bound solutions are possible (see Fig. |3]). 

0.5 1.0 1.5 2.0 2 5 3.0 




Figure 3: 

Bound solutions in three dimensions for r^, = 2. The horizontal line corresponds to k = —5. 

We find the scaling of the total number of bound state solutions with rj, by proceeding 
exactly as in the two dimensional case. Taking the limit \/X — > on the rhs of (fT9l) and 
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using the approximation ( IT2|) one finds that 

i{\f\rb) + K^ ,3 (A/Art) i (Ar, +3) 
hm -2V\^^^ -- = '"^ ^ . (20) 

Equating this to 2k we find that as in the two dimensional case, rim scales linearly with ri, 
in three dimensions also 

.„^-<^. (21) 

The number of bound state solutions (ignoring the degeneracy factor) is then given by 

(2^rfc + 3) 

rib = 1 [rim an mtger) , 

Ub = L'~'^~~^~~~^J (^m ^ot an integer) . (22) 

Taking into account the degeneracy in the spherical harmonics, the total number of bound 
states is given by S*^Q~^(2nj + 1) = . Using the expression for Ub in ( l22l) we see that this 
has terms depending quadratically and linearly on rb as well as a term independent of it. 
However for Vb^ 0, the quadratically dependent term dominates over the other two. 



3 The Nonrotating BTZ Black Hole 

The analysis of the previous section can be done for black hole spacetimes as well, treating 
the black hole horizon as a boundary. One of the main interests would be to see whether 
similar results hold in that case as well. If they do then it would provide a plausible 
explanation for the origin of black hole entropy without requiring a detailed knowledge 
of quantum gravity. The picture would be a realization (at least for a large black hole) 
that the bound states are localized near the horizon and their number (to leading order) 
is proportional to the area of the black hole. 

With these remarks, we now apply the methods of the previous section to the case of 
BTZ black hole [TP]. To simplify calculations we limit ourselves to the case of nonrotating 
BTZ black holes. We want to obtain the Laplacian for a spatial slice in this background. 
The main point of the calculation is to do a self-adjointness analysis of the Laplacian and 
see if it implies mixed boundary condition which, as noted earlier, would introduce a length 
parameter in the problem that can be taken to be of the order of the Planck length. Thus 
we will need to identify the Laplacian operator on the BTZ spacetime, the metric for which 
(in the nonrotating case) is given by 

ds'^ = g^adx^dx^ = - '^^ ~J^^ dt^ + —J-^dr^ + r^dO^ , (23) 

where r_|_ corresponds to the horizon and is the negative inverse of the cosmological 
constant. The massive Klein-Gordon equation is given by 

(□ - m')<P = -^d^i^g'^^dp^ - = , (24) 
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where 

gap ig 

the inverse of the metric in ( 123|) and ^J—g is its determinant. Using (123|) in 

( 12^1) we get 

+ - = . (25) 

This equation has been analyzed in detail at several places, see [20] for instance, and it is 
known that the equation can be solved exactly. 



dU + -dr 



p 



-dr 



3.1 The measure in the r coordinate 

We are interested in the self-adjointness of the Laplacian which requires a suitable measure 
(defining the inner product) on the Hilbert space. In a general curved background the scalar 
product is defined by [21] 

(01,02) = -tJjMd.r,) - (a^0i)0;]y^n^rfS. (26) 

Here is the future directed unit vector orthogonal to spacelike hypersurface S; is 
the determinant of the metric on the hypersurface with dH = drdO being the volume 
element on it. For the BTZ metric, (l23!l . the unit normal to the hypersurface is given 
by g°'^nanp = g^^ntUt = —1. This implies that the only nonzero component of is 
n* = —l/{r'^ — r^)^/^. We also have = Irjir^ — r^)^/^ and ciS = drdd. Thus the inner 
product for nonrotating BTZ is given by 

(01,02) = I [ [Mdt<t>*2) - Wi)0;] , f'' 2, drde. (27) 
Note that the measure diverges near the horizon (this is related to the redshift). 



3.2 Equations in tortoise coordinates 

It turns out that the problem can be translated in the form of standard Schroedinger equa- 
tion by going to the tortoise coordinates. This also makes the analysis of self-adjointness 
issue easier. We use the ansatz = e~*^*e*"'^^^ in Eq. fl25|) and obtain 

(r^-rl)d'^f 2rdf 3/ rlf Pu^f 2 . 



P dr"^ P dr AP ArW - r\ 

To aid further calculations we now go to the tortoise coordinates which are effected by the 
transformation 

P 

P , _i / r 

= coth — 

r+ \r+ 

r = -r+coth(^) . (29) 
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Note that the horizon r = r_|_ in tortoise coordinates is at r,,, = — oo and spatial infinity is 
at r,, = 0. Computing drf and 9^/ and substituting in (128|) we get 



' +17 — + ——-m^M]f + co'f = 0. (30) 



We now transform Eq. (!30l) entirely in terms of n 



f^V 1 /3a2 \ 1 / 2 2^ 

+ -727— T^ + ^'^ / + — 727— T^ + 7T / = ^V, (31) 



(ir2 sinh^(an) V 4 / cosh^(ar^,) V 4 



where we have used the notation a = vM/l (we also note that the horizon is given by 
r_|_ = l\fM, M being the mass of the black hole). This equation is in the form of standard 
Schroedinger equation 



1 /3a2 ^ I (^2 ^2 



with the potential 

V = ^ ( — + m^M ] + ^ ( — + ^ 1 • (32) 

sinh2(«r,) V 4 ; cosh2(ar,) V 4 /V ^ ^ 

We now need to write the scalar product in terms of the r^, coordinate and for this we 
note that in the ansatz (p = e-«'^*e™^'/i^ ^i^^^e is an explicit factor of ^/r sitting in the 
denominator. The two factors of (p in the scalar product will thus contribute a factor 
of which will cancel the r in the numerator of ( !27l) . Thus, in terms of the tortoise 
coordinates the integration measure is simply dr^,. It is now straight forward to verify that 
the Laplacian is a self-adjoint operator for the boundary condition nf = df /dr^ imposed 
at the horizon = —00. This is of the same form as the condition we had in the previous 
section on flat backgrounds. 

Knowing the appropriate boundary condition that makes the Hamiltonian operator self- 
adjoint, we can now go back to the r coordinate in terms of which the boundary condition 
is 

Ag{r) + {r-r+)drg{r) = Q, (33) 

where g{r) is the solution of the radial part of the Laplacian in the original coordinates. 
In the above equation both sides are evaluated at r = r+ and A is a constant depending 
on K, r+, /. 

Having obtained the conditions for self-adjointness in (l33l) the analysis requires numeri- 
cal computations. This is not completed because of strong oscillations of solutions near the 
horizon for generic boundary values. But we will argue that analysis will lead to similar 
results. 

Equation ( l30|) is the Hamiltonian eigenvalue equation in 1-dimension for each value of 
n. Comparing this with similar Eq. (jl]) we flnd that the self-adjointness condition can be 
satisfled only for flnite n up to some maximum value rim ■ The eigenvalues are negative and 
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oc — since contributions come from the second order derivative term in (128!) satisfying 
boundary conditions. These are close to the boundary. From our flat space computations 
it is clear that the number of bound states is finite and proportional to r+ (r^) in two 
dimensions (three dimensions). The same result continues to be true even in the BTZ 
case as can be seen by near horizon analysis of solutions of ( |28l) satisfying the boundary 
conditions f l33|) . though obtaining the exact number of solutions looks difficult. The reason 
for this is that in the Hamiltonian f pOj) . we see that for large r there is a term in the 
potential which is proportional to like the harmonic oscillator. Hence, the spectrum 
has eigenvalues which are positive and wave functions are supported over a larger length 
scale controlled by cosmo logical constant. We should separate the two kinds of states, 
and this can be analyzed in the large r_|_ ^ / ^ 1/k regime. This problem is akin to 
self-adjointness for harmonic oscillator on — D. Furthermore the boundary (horizon), 
which is at r^, = — oo, is actually located at a finite distance r+ in usual coordinates. 
Hence the self-adjointness condition is related to the dilatation operator (r — r^)dr near 
the horizon. The scaling operator along with the Hamiltonian forms elements of an SL(2,R) 
algebra which has been used through conformal quantum mechanics to understand the near 
horizon dynamics [22| |23] . 

4 Discussion 

In this paper we find novel bound state solutions for Laplacian on manifolds with bound- 
ary. These arise by requiring it to be self-adjoint. These states are localized close to 
the boundary and serve to explain the states of a black hole contributing to the entropy. 
They also contribute to Hawking radiation which will be explored later. These features 
arise through brickwall mechanism and quantum mechanical origin justifying the boundary 
conditions on the horizon. The horizon states can be compared with membrane paradigm 
also ES]- They can also be looked as states obtained from spin network of LQG though 
our states contribute directly to energy whereas spin network give area bits. 

In this connection we wish to point out an interesting result from [H]. Here a general 
analysis is done for self-adjoint extensions and the most general boundary conditions are 
characterized by an infinite dimensional unitary matrix linking the boundary data. It is 
given by 



evaluated at the boundary and n stands for the normal. They also point out that if U 
has an eigenvalue —1 then the extensions given by f/e~** have a negative energy state 
which is an edge state. They are characterized by the length parameter. However, when 
the Hamiltonian has in addition potential term which contributes positive energy to the 
state which increases with tt,, then there is a balance of these contributions resulting in 
finite number of states. This results in the number of edge states being constrained by the 
radius. 

In the case of Schwarzschild black hole the situation is different since it is asymptotically 
fiat. This work is in progress and will be presented elsewhere. 




(34) 
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